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Abstract 

We study linear difference equations with variable coefficients in a ring using a new non- 
linear method. In a ring with identity, if the homogeneous part of the linear equation has a 
solution in the unit group of the ring (i.e., a unitary solution) then we show that the equation 
decomposes into two linear equations of lower orders. This decomposition, known as a semi- 
conjugate factorization in the nonlinear theory, generalizes the classical operator factorization 
in the linear context. Sequences of ratios of consecutive terms of a unitary solution are used 
to obtain the semiconjugate factorization. Such sequences, known as eigensequences are well- 
suited to variable coefficients; for instance, they provide a natural context for the expression 
of the classical Poincare-Perron Theorem. We discuss some applications to linear difference 
equations with periodic coefficients and also derive formulas for the general solutions of linear 
functional recurrences satisfied by the classical special functions such as the modified Bessel and 
Chebyshev. 



1 Introduction 

A linear, non-homogeneous difference equation with variable coefficients is defined as 

X n +1 = aO,nX n + Ol jn X„_l H h CL k , n Xn-k + K (1) 

where {b n } and {dj,n} are given sequences in a nontrivial ring R for j = 0, 1, . . . , k. 

Equation (pQ) may be unfolded in the standard way to a map of a (k + l)-dimensional i?-module 
(or vector space, if R is a field). In this form there is substantial published research that extends 
the classical theory on the real line to a diverse selection of rings that include finite rings, rings of 
polynomials and other rings of functions as well as abstract rings of various types. See, e.g., [1], 

0, 0, m, m, ph, eel \m- 

We pursue a different approach in this paper, using a method originally developed for nonlinear 
difference equations. Specifically, we explore the existence of a semiconjugate factorization of ([I]) 
in its underlying ring R into two or more linear difference equations of lower orders. While a linear 
equation may have many different semiconjugate factorizations, we seek one where the equations of 
lower orders are also linear. The key requirement for the existence of a semiconjugate factorization 
is the existence of a special type of sequence that we call an eigensequence of the homogeneous part 
of ([!]); eigenvalues are essentially constant eigensequences. 
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We show that in a ring R with identity, quotients of the consecutive terms of a unitary (invert- 
ible) solution form an eigensequence in R. So it is not surprising that the classical Poincare- Perron 
Theorem finds a simple and natural expression in this context; see Section [7] below. Semiconjugate 
factorization generalizes the classical notion of operator factorization of a homogeneous equation 
with constant coefficients over the real or complex numbers. More generally, semiconjugate factor- 
izations of linear equations over arbitrary nontrivial fields have been studied in Chapter 7 of |19| . 
The results in this paper substantially extend those in [19] to rings and add several new results, 
including results for equations with periodic coefficients and for recurrences in rings of functions; 
see Sections l5lfT0l below. 

Establishing the existence of an eigensequence and calculating it are problems that are more 
challenging in rings than in fields but the rewards are also potentially greater. In particular, rings of 
functions on a given set typically are not fields but the semiconjugate factorization method applies 
to linear difference equations in such rings in essentially the same way (finding an eigensequence) 
as it does to linear difference equations in a field. We study semiconjugate factorizations of linear 
difference equations in rings of functions and apply the results to obtain explicit formulas for the 
solutions of some known higher order functional difference equations; e.g., second-order recurrences 
that define modified Bessel functions or Chebyshev polynomials, but with arbitrary initial functions. 

2 Preliminaries 

A (forward) solution of (pQ) is defined, as usual, to be any sequence {x n } in R that satisfies the 
equation for n > 0. Given the recursive nature of ([1]) it is clear that with any k + 1 given initial 
values Xj £ R for j = 0, 1, . . . , k (the number k + 1 being the order of the difference equation) 
(TjQ) generates a unique solution in R through iteration, since R is closed under its addition and 
multiplication. If b n = for all n then ([T]) is homogeneous and in this case, the constant sequence 
x n = for n > is a solution of ([I]), namely, the trivial solution. 

In the classical theory of higher order linear difference equations in the field of real numbers, 
linear operators such as the ones used to define ([TJ may be "factored" using their eigenvalues; see, 
e.g., Section 2.3 in [5]. This elementary procedure yields both a reduction of order for a linear 
difference equation and a symbolic "operator method" for obtaining its solutions. For discussions 
of these basic classical notions, including operator methods, see [5] or [9]. 

In this section, for the reader's convenience we present some general results from |19j that are 
valid for all difference equations of recursive type, not just the linear ones. 

Let Q be a nontrivial group and consider the recursive difference equation, or recurrence 



where f n : Q + — > Q is a given function for each n > 0. Starting from a set of k + 1 initial values 
Xj £ Q for j = 0, 1, . . . , k a unique solution of (J2]) is obtained by iteration. 




(2) 



2 



Equation ([2]) may be unfolded in the usual way to a first-order recurrence 

Xn+i = $n(X n ) 

on Q k+1 where $ n : Q k+1 — > Q k+1 . Let k > 1, 1 < m < k. Suppose that there is a sequence of maps 
: G m ~^ G m an d a sequence of surjective maps H n : Q k+1 — > Q m that satisfy the semiconjugate 

relation H n+ \ o J„ = o H n for a given pair of function sequences {3n} and {& n }. Then we say 

that 3« is semiconjugate to for each n and that the sequence {H n } is a /orm symmetry of ([2]). 

Since m < A; + 1, the form symmetry {H n } is order-reducing. 
We state the next core result from [19] as a lemma here. 

Lemma 1 (A semiconjugate factorization theorem) Let Q be a nontrivial group and let k > 1, 

1 < m < k be integers. If h n : g k - m + 1 — ^ Q i$ a sequence of functions and the functions H n : 
gk+i _^ gm are define fry 

H n (u , ... ,u k ) = [u * h n (ui, . . . 

where * denotes the group operation in Q then the following statements are true: 

(a) The function H n is surjective for every n > 0. 

(b) If {Hn} is an order-reducing form symmetry then the difference equation (dP is equivalent 
to the system of equations 

tn+l = tpnitm ■ ■ ■ > tn— m+l)i (3) 
Xn+1 = tn+l * K+l(x n , ■■ ■ , Xn-k+rnY 1 (4) 

whose orders m and k + 1 — m respectively, add up to the order of 

(c) The map <3> n : Q m —> Q m is the standard unfolding of Eq.^) for each n > 0. 

Remark 2 Part (c) above permits us to stay within the context of higher order difference equations. 
Being able to work within this context is especially beneficial in the case of nonrecursive equations 
(including some linear ones) which do not in general unfold to maps on modules over rings (or 
vector spaces over fields) and therefore, their solutions are not determined via group actions; see 
the final section of this paper for an illustrative example. Extensions of the method of this paper 
to nonrecursive equations no longer rely on semiconjugacy but they do retain the basic concepts of 
form symmetry and factor- co factor pairs; see U9\j . Chapter 8. 

Definition 3 The pair of equations (Qp and ^ constitute the semiconjugate factorization, 
or sc- factorization of (Hp. This pair of equations is a triangular system (see [22]) since is 
independent of |^p. We call ^ the factor equation of and its co factor equation. 



3 



Note that (|3|) has order m and has order k + 1 — m. Consider the following special case of 
H n in Lemma Q] with m = k 

H n (u Q ,ui, . . . ,Uk) = [uq * h n (ui),ui * h n -i(u 2 ), . . . ,u k -i * h n -k+i(uk)\ (5) 

where h n : Q —> Q being a given sequence of maps. The semiconjugate factorization of (j2J) in this 
case is 

tn+l = 4>n{tn, ■ ■ ■ , £n-k+l), (6) 

Xn+1 = t n +l * K+liXn^ 1 (7) 

in which the factor equation has order k and the cofactor equation has order 1. 

The next result gives a necessary and sufficient condition for the existence of a form symmetry 
of type (|5J); see [19] for the proof. 

Lemma 4 (Invertible-map criterion) Let Q be a nontrivial group and assume that h n : Q — > Q 
is a sequence of bisections. For arbitrary elements Uq, vi, ■ ■ ■ , G Q and every n > define 
(o,n(uo) = uo and for j = 1, . . . ,k define 

(j,n(u ,Vi,. . .,Vj) = h~l j+l ((j-l jn (u ,Vi,.. . *Vj). (8) 

Then {IjJ has the form symmetry |5jj if and only if the quantity 

fn((0,n, Cl,n(^0, Vi), . . . , Cfc,n(«0, "1, • • • , «&)) * K +1 (u ) (9) 

is independent of uq for every n > 0. In this case has a semiconjugate factorization into 
(0|) and ([?[) where the factor functions in (0|) are given by 

</> n (vx,...,V k ) = /n(Co,n,Cl,n(^0,Wi),...,Cfe,„(MO,«l,---,Wfc)) * h n+1 (u ). (10) 

In the context of rings, the group Q in the preceding result is the additive group of the ring so 
that * denotes addition and thus (jSj), © and ([7|) read, respectively, as follows 

Cj,n(U0,Vl,-. ■ ,Vj) = h~^ j+1 (Vj - Cj-l,n( u 0,Vl, ■ ■ -,Vj-l)), 
fn((o,n,(l,n(u Q ,Vi), . . . ,(k,n(lM),Vi, ■ ■ ■ ,«fc)) + K+l(uo), and: 

Xn+l = — hn+1 [x n ). 

A basic class of maps h n in rings is defined next. 

Definition 5 Let R be a ring and let {a n } be a sequence in R such that a n ^ for all n. A linear 
form symmetry is defined as the special case of ([5|) with h n (u) = —a n u; i.e., 

[uq - a n ui,ui - a n -iu 2 , Uk-i - a n ^ k+1 u k ] (11) 



4 



If a is not a zero divisor then h{u) = —au is one-to-one or injective since for every u,v £ R 

h(u) = h(v) => — a (it — u) = 0=^u — u = 0. 

In general, h is not surjective even if R contains no zero divisors (consider a £ Z, a / But if 
R has an identity and each a is a unit then each h is a bijection with inverse h~ l (u) = —a~ l u. 

Remark 6 Many nonlinear difference equations possess the linear form symmetry $11]) : see J18\j 
and \Wj. 

3 SC-factorization in rings 

Assume that the underlying ring R of (pQ) has a (multiplicative) identity denoted by 1. For such a 
ring, the set of all units (elements having multiplicative inverses, or reciprocals) is a group, namely 
the unit group, that we denote by G. 

A unitary sequence is any sequence in G. If {u n } is a unitary sequence then the sequence 
{un+iu^ 1 } of right ratios of {u n } is well-defined and unitary. Similarly, the sequence {u~ l u n +\} 
of left ratios is well-defined and unitary. If R is commutative then the sequences {u n+ \u~ 1 } and 
{u~ l u n+ \\ are the same, representing the ratios sequence of {u n }. 

Call two sequences {x n } and {y n } in R right equivalent (or left equivalent) if there is a unit u 
such that y n = (or y n = ux n ) for all n. These two relations on the set of sequences in R are 
indeed equivalence relations, and if {x n } is unitary then so are {x n u} and {ux n }. 

The next result has the same flavor as the result in calculus which states that differentiable 
functions having the same derivative are equal up to a constant. 

Lemma 7 Let R have an identity and {x n } and {y n } be unitary sequences. Then {x n } and {y n } 
are right (or left) equivalent if and only if their sequences of right (or left) ratios are equal. 

Proof. Suppose that {x n } and {y n } are right equivalent. Then y n = x n u for some unit u and all 
n so that 

Vn+lVn 1 = X n+ iu(x n u) _1 = X n+ iu{u^ X^ 1 ) = X n+ iX~\ 

i.e., the sequences of right ratios are the same. Conversely, suppose that T/n+iy" 1 = Xn+ix^ 1 for 
all n > and define u = Xq 1 yo. Then 

xiu = xiXQ l y = yiy^yo = y\. 

This equality also implies that u = x^ X y\ so the preceding argument may be repeated to show 
that x n u = y n for all n > 0. A similar argument proves the left-handed case. ■ 

Recall that the existence of a linear form symmetry (llip implies that ([2]) has a semiconjugate 
factorization with a first-order, linear non-homogeneous cofactor equation 

+ a n+1 x n . (12) 
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The following necessary and sufficient condition for the existence of a linear form symmetry is 
a consequence of the invertible-map criterion with Q being the additive group of the ring. 

Lemma 8 Let R be a ring with identity. Equation has the linear form symmetry if and 
only if there is a unitary sequence {a n } in R such that the quantity 

fnim, Cl,n(uO, Vl), . . . , (k,n(uQ, V\, . . . , V k )) ~ «n+1^0 (13) 

is independent of uq for all n, where for j = 1, . . . , k, 

Q >n {uQ,vx,. . . ,vj) = I JJa^-j J n -^ JJa n _ m+ i «i. (14) 
\i=0 / i=l \m=i J 

Proof. Define h n (u) = —a n u for each n so that h~ l {u) = —a~ l u for all n. If we define Cj,o = 
and for j = 1, . . . , k set 

0,n(«0,Vl, -->,Vj) = a~^ j+1 [Cj-l,n(uo,Vl, ■ ■ -,Vj-l) - Vj]. 

recursively then the first assertion of the lemma is true by the invertible-map criterion. To prove 
(f!4"|) . observe that 

Ci,n(wo, vi) = a~ l (u Q - vi) = a~ x u§ - a~ V 



which proves fjl4[) if j = 1. Suppose that (|14p is true for indices less than j where j < k. Then 



Cj,n(«0,«l, ■■■,Vj) = a n ^ j+1 [Cj-l,n(uO,Vl, ■ ■ ■ ,Vj-l) ~ Vj] 



a 



-l 

n-j+l 



(3-1 

IT oin - i 

\i=0 

/i-i 

n ° in - i 

Vi=0 



'i-2 

vl=0 

-1 



3-1 1 3-1 N 
^0 - JJ"n-m+l 



^ - 



i=l \m=i 



-1 



-1 



3-1 / 3 

U ° ~ Yl [ II a n-m+l I «i - <*n-j+l v 3 
i=l \m=i / 

1 A , A X "I 



WO - ^ On-m+1 



i=l \m=i 



and the proof is complete. ■ 

The following application of Lemma [8] gives a sc-factorization theorem for linear difference 
equations in a ring with identity; also see the comments following the theorem. 
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Theorem 9 Let R be a ring with identity. The linear equation (OP has the linear form symmetry 
Ul\) with unit coefficients if there is a unitary sequence {a n } that satisfies the relation 

k \ _1 

Oi n +l = «0,n + ^2 CLj : n I Oin-i J . (15) 
j=l \i=o J 

The corresponding semiconjugate factorization of ([!]) is 

H 1" 4-l,^n-fc+l + bn (16) 

Xn+i = a n+ ix n + t n+ i (17) 



where for m = 0, . . . , k — \, t m+ \ = x m+ \ — a m+ ix m and 



-l 

k I i \ 

J 



®"in,n 



i=m+l \ j=m+l 



Proof. By Lemma [8] it is only necessary to determine a unitary sequence {a n } in R such that for 
each n (|13|) is independent of u$ for the following functions 

f n (u , ...,u k ) = a 0tn u + a 1>n ui H h afc^fc + K- 

For arbitrary uo, u i, • • • , v k £ R and j = 1, . . . , k define Cj,n(uo, v\, . . . , Vj) as in LemmaO Then 
expression (fT5|) is 



-a n+ lU0 + &n + a0,n^0 + ai t nCl,n(uo,Vi) H h ak,n(k,n{u , Vl, . . . , Ufc) 



6n + 



/i-i 



-Q 



n+1 + °0,n + ^ aj,n I ]^[ 



Or 



i=0 



k 3/3 
UQ - %n ^ On-m+1 



i=l \m=i 



The above quantity is independent of uo if and only if the coefficient of uo is zero for all n; i.e., if 
and only if {a n } satisfies the difference equation (|15|) , Dropping the uq terms leaves the following: 



bn ^ ®j t n 

3=1 



3/3 

EE 

8=1 \m=j 



-1 



n 



n— m+l 



-r 



i=j 



(18) 



From this expression, we obtain the sc- factorization of J!]). The cofactor equation (|17j) is simply 
(|12p while the factor equation is obtained using ©, the above calculations and (|18p . Finally, (|16p 
is obtained by slightly adjusting the summation indices to simplify notation. ■ 
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Remark 10 (The inversion form symmetry) Equation U5\) is not only a consequence of the in- 
vertible map criterion but it is also related to a different semiconjugate factorization. Think of{a n } 
as a solution of the following k-th order difference equation on the (multiplicative) unit group G of 
R 

II r »-* • (19) 

j=l \i=0 J 

Consider the homogeneous part of |7]) i.e., the linear equation 

x n +i = ao, n x n + ai n x n -i H h a k , n x n -k- (20) 



If a j,n £ G for all j then $19\) turns out to be a factor equation of a sc- factorization in G of i20\) 
and the cofactor equation is x n+ \ = r n+ \x n ; see }2(Jf . The inversion form symmetry that yields this 
sc-factorization is [uqU^ , U1U2 , ■ ■ ■ , Uk-iu^ ]. This (nonlinear) form symmetry is characteristic 
of all difference equations that are homogeneous of order one, linear or not; see Chapter 4 in- 119]. 
Notice that the factor equation U9\) of h20\) relative to the inversion form symmetry is nonlinear 
whereas the factor equation [Tb)) relative to the linear form symmetry is linear. This makes a further 
reduction of order more difficult for l\19\) which is also not homogeneous of order one. 



4 Characteristic equation and eigensequences 

The sequence {a n } in Theorem [U] is a solution of the nonlinear difference equation (|19p . Since {a n } 
plays a fundamental role in the sc-factorization of JT]), it is necessary to examine (|19p closely. This 
equation may be written in a way that does not involve inversion. Multiply it on both sides by the 
quantity r n r n _i . . . r n _ k+l 



r n+lfn r n—l • • • f"n— k+1 — O0,n 

{r n r n _i . . . r n _ fe+ i) + oi )Tl (r n _i . . . r n ^ k+1 ) + 

+ 0,2,n(r n -2 ■ ■ ■ r n „ fc+1 ) + • • • + a k -l t n r n-k+l + flfc,n 

which may be written more succinctly as 

k k-l (k-\ \ 

Y\ r n - i+ i - ^2 a j,n Y\ Tn ~ i ~ afc - n = °- ( 21 ) 

i=0 j=0 \i=j J 

This equation is not as esoteric as it may appear at first glance. To clarify, consider the special 
homogeneous case with constant coefficients, i.e., 

x n+ i = a x n + aix n -i H h a k x n _ k . (22) 
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Then (|2ip reduces to the following difference equation 

r n+1 r n . . . r n _ fc+ i - a (r n . . . r n _ fc+ i) - ai(r n _ 2 . . . r n -k+l) dk-lTn-k+l - a k = 0. (23) 

A constant solution (or fixed point) r n = r of (|23p must satisfy the polynomial equation 

r k+1 - a r k - a ir k ~ l a k = 0. (24) 

The right hand side of (|24p is recognizable as the characteristic polynomial of (|22 [) whose roots 
are indeed the eigenvalues of the linear homogeneous equation ([22]) . 

Definition 11 T/ie difference equation \21}) in a ring R is the characteristic equation of the 

homogeneous part of f7]), i.e., the linear difference equation A20\). Each solution of h21\) in R 
is an eigensequence of A20\) . An eigenvalue is a constant eigensequence. An eigensequence 
whose every term is a unit in the ring is unitary. An eigensequence containing a zero divisor is 
improper. 

Equations ([2T]) and (fl~9j) are not equivalent: every solution of (fl~9j) is a unitary solution of ([2T]) 
but nonunitary solutions may also exist for (|2ip that translate into eigensequences that are not 
unitary, or not proper. 

Example 12 The second- order linear difference equation 

x n+ i = 2x n - 4x n _i, x , xi € Z (25) 

has the characteristic equation r n+ ir n — 2r n -\-A = 0. This has no constant solutions (eigenvalues) in 
Z since the polynomial r 2 — 2r + 4 has complex roots r = 1 ± i\/3. But this characteristic equation 
does have a (nonunitary) period 3 solution given by {1,-2,4,1,-2,4,...} as may be checked by 
direct substitution. This eigensequence is trivially proper since Z has no zero divisors and it is in 
fact unitary in the field Q of rational numbers. 

The difference equation \25\) is also valid in finite rings Z m of integers modulo a given positive 
integer m and different cases occur. For instance, TLyj is a field so the above eigensequence is 
unitary. But in 7L\% the same sequence of period 3 is an improper eigensequence since all even 
numbers are zero divisors. Improper eigensequences are not desirable because ring extenstions do 
not render them unitary and Theorem cannot be applied. Finally, we mention that for some 
values of m the polynomial r 2 — 2r + 4 has roots in Z m which are eigenvalues of &25\); e.g., in the 
field Z5 the roots are 3 and 6 while in the ring Z12 the roots are —2 = 10(modl2) and 4 both of 
which are improper. 

Example 13 Consider the difference equation 

x n +i =x n + X n _i (26) 
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that is also known as the Fibonacci recurrence because with initial values xq = and x% = 1 Q2t 
generates the Fibonacci sequence 1,1,2,3,5,8,13. .. , denoted {F n }. The characteristic equation of 
is 

r n +ir n ~ r n - 1 = 0. (27) 
This equation has no solutions in the ring of integers Z ; constant or otherwise. For let ri, ri € Z 



and note that n 7^ because clealry r n = (foes no£ so/we |i?7| ), iVou; (£7) /ias a solution r^ £ Z i/ 
and on/y i/ri = ±1. Either T\ = 1, r2 = 2 so i/iai r3 = 3/2 Z or ri = —1, r2 = and no wa/ne 
is defined for r%. Hence, 126\) has no eigensequences in Z. 

T/ie eigenvalues (constant eigensequences) of this equation are roots (1 ± \/5)/2 0/ i£s charac- 
teristic polynomial r 2 — r — 1. Thus \20) has no eigenvalues in Q; but unlike Z, in Q fj?7| ) can be 
stated as r n+ i = 1 + l/r n . Iteration starting from (say) rg = 1 yields r n = F n+ i/F n , a unitary 
eigensequence for \26) . Theorem then yields a sc-factorization, in Q, of $26\) consisting of the 
pair of equations 

tn+l = a o,n^rc; a 0,n 



F n 



+1 

_ Fn+2 , 
Xn+1 — j-, Xn 1 In+1- 

We close this section with the following corollary of Theorem [9] for second order equations. See 
the next section for examples. 

Corollary 14 Let R be a ring with identity and assume that the difference equation 

has a unitary eigensequence {r n } in R; i.e., r n+1 = a 0n + a\ n r~ l for all n > 1. If there is an 
integer m > 1 such that a\^ m = then 

x n = r„r n _i • • • r m+1 x m for all n > m + 1. 

Proof. By Theorem [9] the second-order difference equation has a sc-factorization 

tn+l = a 0,n^n = ~ a >l,n r n ^ni t\ = X\ — r±XQ 
Xn+1 — fn+\.Xfi -\- t n +l 

If ai >m = then t m+ \ = 0. Thus, t m+ 2 = -am+i^m+iim+i = and by induction, t n = for 
n > m + 1. This implies that the cofactor equation reduces to x n = r n x n -\ for n > m + 1. Upon 
iteration we obtain the formula for the solution {x n } of the second-order equation. ■ 
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5 Periodic coefficients 

In this section we study the following difference equation with periodic coefficients in a nontrivial 
ring R, i.e., 

■Eji+l — Q"n&n ~^ bnX n —i, Q, n + Pl — (2 n , & n -|_p 2 — b n , n — 0, 1, 2, . . . (28) 

where the (minimal or prime) periods pi,P2 are positive integers with least common multiple 
p = lcm(pi,p 2 ); a central question is whether (|28p has an eigensequence of period p in R. If 
so then such an eigensequence yields a sc-factorization of the second-order equation into a pair 
of first-order ones. This may occur subject only to algebraic restrictions on the ring concerning 
the existence of roots for a quadratic polynomial but regardless of whether (|28p has any periodic 
solutions. In general, one cannot expect an equation such as (f28|) to possess any periodic solutions 
but the problem of existence of peridoic solutions has been studied previously; see e.g., [H] . 

An eigensequence of period p exists in R if there is an initial value r% 6 R such that the 
characteristic equation of (|28l) . i.e., the first-order quadratic difference equation 

r n+1 r n = a n r n + b n (29) 

has a solution of period p in the ring R. Suppose that there are rj € R that satisfy (|29p for 
j = 1,2, ... ,p. Then 

r 2 ri = a\r\ + b\, r 3 r 2 = a 2 r 2 + b 2 

Let L\ = r\ so that r 2 L\ = a\L\ + b\. For j = 2, . . . ,p define Lj + \ = ajLj + bjLj^i. Then 
r 2 r\ = r 2 L\ = L 2 so that 

?"3^2 = (7~3r 2 )ri = a 2 r 2 ri + b 2 n = a 2 L 2 + b 2 L x = L 3 , 
7*4^3 = {rAr 3 )r 2 ri = a 3 r 3 r 2 r 1 + b 3 r 2 ri = a 3 L 3 + b 3 L 2 = L 4 , 

By induction, for j = 2, . . . ,p 

r i+i L J = ( r i+i r i) r i-i ■■■ri = ajrjTj-i • • • n + 63^-1 ■ ■ ■ n = a^Lj + fc^-i = L i+ i. 
This process yields a solution {r n } of (f29|) with period p if and only if r p+ i = n; thus, 

riLp = r p+ iL p = a p L p + b p L p _i (r x - a p )L p = 6 p L p _i (30) 

The quantities L\,...,L p that are generated above evidently depend on r\ in a linear way so 
there are a.j , f3j £ R such that 

Lj = a j r 1 + (3j 
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for j = 1,2, ... ,p. Inserting this form in (|30p yields 

(n - a p )(a p ri + (3 P ) - b v {a v -\r\ + /3 p _i) = 
ria p ri + ri/? p - (a p a p + 6 p a p _i)ri - (a p /3 p + 6 p /3 p _i) = (31) 

The definition of Lj implies 

aj+iri + /3 j+ i = aj(ajn + fy) + bj(aj-\r\ + 

= (a/ay + bjaj-i)ri + o^-jSj + bj/3j-i 

Suppose that R has an identity 1. By matching coefficients on the two sides of the above 
equality, we see that the coefficients otj,/3j satisfy (|28p for j = 1, 2, . . . ,p with initial values 

a = 0, ai = 1; /3 = 1, jSi = 0. (32) 

Using this fact to simplify (|31 1) we conclude that if n is a root of the following polynomial 

ra p r + r/3 p - a p+ ir - fip+i = (33) 

then the solution {r n } of (|29p has period p. These observations prove the following result. 

Theorem 15 Let R be a ring with identity 1 and for j = 1,2,... ,p, let ctj,/3j be obtained by 
iteration from [28\) subject to [32\) . 

(a) If r\ is a root of the quadratic polynomial (E3|) in R and there are rj £ R satisfying l2§\) for 
j = 2,... ,p then {r n }^ =1 is an eigensequence of [28\) with preiod p. 

(b) If a root r\ of $33\) in R is a unit and the recurrence 

r j+1 = aj + bjrj 1 (34) 

generates units r-i-, ■■■ ,r p in R then {r n }^ = i is a unitary eigensequence of A28\) with preiod p that 
yields the sc- factorization 

t n +i = -b n r~ l t n , ti= x 1 - nx 

The polynomial (|33p simplifies further if the coefficients aj,bj are in the center of R. Then 
a.j , f3j are also in the center of R so (|33l) reduces to 

a p r 2 + (f3 p - a p+ i)r - (3 P+1 = 0. (35) 

Note that if a n = a and b n = b are constants then (with p = 1) the quadratic polynomial (|35p 
takes the form r 2 — ar — b = 0, i.e., the characteristic polynomial of the autonomous linear equation 
of order 2. 
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Example 16 Consider the difference equation 



x n _i_i = 2 cos 



2-irn 



x n + X n —\ 



(36) 



where a n = cos27rn/3 has period 3 with a\ = a 2 = —1 and a% = 2 and b n = 1 is constant (period 



1). Let us assume that the underlying ring is the field 
readily calculated from [36]) using HM): 



a>2 



of real numbers. The numbers cxj,/3j are 



-1, a 3 = 2, a 4 = 3, (3 2 = 1, h 



-1, /3 4 



The quadratic equation t35\) with p = 3 is 2r 2 — 4r + 1 = in this case. Of the two real 
roots (2 ± y2)/2, let r\ = (2 — \/2)/2. T/ien using {3$ we readily calculate r 2 = 1 + v2 7 r 3 = 
—2 + \/2. Since these are units in R, 6y Theorem [75] an eigensequence with period 3 is obtained. If 
p = —l/(r\r-2Tz) then the sc- factorization of {36}) is obtained by straightfoward calculation with the 
factor equation 



%+i — p'h, t 3 j 



+2 



phi 
n 



> *3?+3 



rir 2 : 



i > 0, *i = X\ - T\Xq 



and the cof actor x n +i = r n+ \x n + t n +\. Since p = 1 + \/2 aW solutions of \36\) with t\ ^ are 
unbounded. However, t\ = w/ien */ie initial values satisfy x\ = r±xo- Then t n = /or a// n and 



X3r, 



(-l) n X 



£3n+l 



ri(-l) n x 



^3n+2 



rir 2 (-l) n x 



, n > 1. 



These special solutions of \36\) converge to exponentially for all xq. 

The above calculations are meaningful in rings other than R. By way of comparison, now 
suppose that the underlying ring of \30) is a finite field of type Z p where p > 5 is a prime. In 
particular, the quadratic polynomial 2r 2 — 4r + 1 factors in TL-j with roots 3, 6. If n = 3 then using 
{3$ we find that r 2 = 4 and r$ = 1 do in fact yield an eigensequence with period 3. Then, modulo 
7, p = — l/(rir 2 r3) = 4 and 



„3i 



V 



1, P 



3i+l 



n 3i+2 



and */ns pattern of period 3 for distinct powers of p yield a pattern of period 9 for t n as shown in 
the following table 



i 





1 


2 


hj+i 


*i 


4ti 


2ti 




2ti 


*i 


4*i 


t3j+3 


3ti 


5ti 


6*i 



where t\ = x\ — 3xo- /n calculating the above entries we used l/r\ = 5 and l/rir 2 = 3 (modulo 7). 
Of course all solutions in Z7 aZZ {*„} and {x n } will be periodic with period at most 48 since Z7X 
Z7 has that many points (besides (0,0) which yields the trivial solution). 
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If quadratic polynomial (|33p (or (|35p in the commutative case) has no roots in the underlying 
ring then periodic eigensequences with peirod p do not exist in that ring. However, other periodic 
eigensequences may exist; e.g., for the autonomous equation (f2"5"j) where p = 1 recall that there are 
no integer eigenvalues but there is an eigensequence of period 3 in Z. In some cases, a nonperiodic 
but still useful eigensequence may exist as in the next example. 

Example 17 Consider the following variant of A36\) in M: 

x n+ i = 2cos ^-y-^ x n ~ %n-i- (37) 

We find that a<i = —1, a 3 = 0, 04 = 1, 02 = —1, @3 = 1, Pa = 3. With these coefficients 
[So]) has no roots so (3~7\ ) has no period-3 eigensequences. But the recurrence P^] ) can be used to 
generate other types of eigensequences. For instance, if r\ = 1 then it can be verified by induction 
that 

q • 1 1 3j + 2 1 . . n 

rsj+i = 3j + 1, r 3j+2 = — ^j-j-y, r 3j+3 = -gT^' 3 > 
is a (nonperiodic) eigensequence for (37\ ). Note that r n r n+ ir n+ 2 = 1 for all n so 

hj+i = h TT — = h TT = ii = x\ - x 

and the factor equation may be expressed as 

h h 
hj+i - h, h j+2 - - - , t 3j+3 - - — , j > 0. 

T/ie cofactor can now be specified as follows 

x-ij+i = (3j + l)x 3 j 

3j + 2 *! . . 

^3j+2 = - - 1 ^3j+i + ~ 1 = ~(3i + 2)x 3j - t\, 

^+3 = -3^2^+2 - 3^ = x 3j 

The last equation implies that x^j = xq for all j so the general solution of {31\) is 



X-, 



x ifn = 3j 

uxq +xi if n = 3j + 1 

— nxo — x\ if n = 3j + 2 



In particular, if xq = i/ien i/ie solution {0, xi, — xi, . . .} 0/ p?7|) has period 3 for all x% 7^ 0. 
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6 Unitary solutions and eigensequences 

A potential difficulty in applying Theorem [9] is finding the sequence {a n } that satisfies (|19|) . i.e., 
finding a solution of (|21|) . Even if such a solution exists then in many cases finding it directly from 
(f2Tj) is usually not easy. Fortunately, it is often possible to calculate {a n } indirectly, by extracting 
it from a unitary solution of (|2L)p in R; i.e., a solution of (|2(jp that is contained in the unit group G. 
Let {x n } be such a unitary solution for a given set of initial values Xo, x\, ■ ■ ■ ,Xk G G. Multiplying 
(|20p by x~ l and rearranging terms gives 

Xn+\X n — OjQ,n ~\~ ^1 .n'-^n- -1 X n ^2,nXn—2X n ' ' ' ^k,nXn—kX n 
= a ,n + ai.nXn^iX^ 1 + a2 t nX n -2{x^\x n -i)x^ 1 H 

~\~ &k,n%n—k{% n _} ; j r \%n~k+l^(% ri ^} t j r 2'^ ri ~k+2) ' ' ' (x n _ \X n — l)X n 
= a ,n + ai^XnX^Y 1 + a 2t n{Xn-lXn-2Y 1 {XnX^xY 1 H 

+ a k,n( x n— k+l x n — fc) { x n— k+2 x n _j t j l _-\) ' ' ' (^n-^n— l) 

If r n = inS'lj for each n then the above equation can be written as 

r n+ i = a , n + ai^r" 1 + a 2 ,„r~l 1 r- 1 + • • • + a &in r~l fe+1 r^ fe+2 • • • r~Vn X , or: 

r n+ i = a 0i „ + ai^r" 1 + a 2 ,„(r n r n _i) _1 H h a k , n {r n r n _i . . . r n „ A . +1 ) _1 (38) 

which is precisely Equation (fT9l) . Thus, the sequence {r n } of right ratios of {x n } satisfies (flUj) . It 
is often easier to find a unitary solution of (|20p than to look for a particular solution of (|19p . Once 
a unitary solution of (f20l) is identified, an eigensequence may be extracted from it using the next 
result that supplements and completes Theorem [9] 

Theorem 18 Let R be a ring with identity. A (unitary) sequence in R is an eigensequence of i20\) 
if and only if it is the right ratio sequence of a unitary solution of Ii20\) . 

Proof. Let {r n } be a unitary eigensequence of ([20]) . choose xq S G and define Xj = TjXj-\ for 
j = 1, . . . , k. Then Xj G G for each j and 

r j+ iXj = (a , n + ai^r" 1 + a 2 , n rjl 1 rj 1 H h ak, n rj\ +1 rj\ +2 ■ ■ ■ rj^rj^xj 

= O-0,nXj + ai t n r j Xj + 02,n r j — i r j Xj + • ■ ■ + ^fc,"^'— fc+lTj— ' ' ' r j — l r j X j 

= a , n Xj + ai^Xj-t + a2 i „r7 1 1 x i _i H h afc, n ^7-JHl r 7-fc+2 ' ' ' r ]-2 r ]-i x j-^ 

where we used the fact that r^Xj = Xj-\. Similarly, rJ_^Xj-i = Xj-2 which yields a further 
reduction 

r j+ iXj = a , n Xj + ai >n Xj-i + a 2 , n Zj-2 H h a M r 7-fc+i r 7-fc+2 " " " r 7-2 x i-2- 



15 



Next, rj_ 2 x j-2 = x j-3 an d the above calculation may be continued to ultimately yield 

rj+iXj = a 0>n Xj + ai )n Xj-i + a 2 , n Xj-2 H h a k ^ n Xj^ k 

Define the right hand side as Xi+i, then proceed to rj+2&j+i an d repeat the calculate to generate 
a new value 

Xj+2 = r j+2 Xj + i = a , n Xj+i + ai. n Xj + a 2 , n Xj-i H h a kjn Xj +1 ^. k 

The values of x n generated by the above construction satisfy the linear equation (|20|) for n = 
j + 1, j +2, . . . Therefore, {x n } is a unitary solution of (|2"0j) whose right ratio sequence is {r n } (by 
construction). The converse is true by the definition of eigensequence and the argument preceding 
this theorem. ■ 

For the homogeneous difference equation (|22p with constant coefficients the following is true. 

Corollary 19 Let {a{\, i = 1, . . . ,k be constants in a ring R with identity such that a k ^ and 
let G be the unit group of R. 

(a) Equation h22]) has a linear form symmetry if and only if it has a unitary solution {u n } in 
R. In this case, the corresponding sc- factorization is determined by the eigensequence {a n } of right 
ratios a n = u n u^_Y of {u n } as follows 

t n +i = a'otn + «i*n-i + • • • + a' k _it n -k + i, (39) 
where a' m = - J2i=m+i a * (ll}= m +i On-j+l) for m = 0, . . . ,k - 1. 

(b) Every fixed point of \23\) is an eigenvalue of 122]) in R. Such a fixed point exists if and only 
if the characteristic polynomial r k+1 — a$r k — ■ ■ ■ — a k _\r — a k has a root in R. 

Remark 20 Equation ( fggl) is once again linear but with order less than \22]) by one. If also 
possesses a unitary solution in R then it is reducible in order by a second application of the preceding 
corollary. This process may be repeated k times to yield a triangular system of k + 1 first order 
linear equations; see [19]. In particular, if R is an algebraically closed field then the characteristic 
polynomial of [22]) factors completely in R and thus, a set of k + 1 eigenvalues is made available 
for a full sc-factorization. 

Example 21 Consider the following third-order linear difference equation in the finite field Z p 
where p is a prime number 

x n+ i = 2x n -i + x n -2, x , xi,x 2 £ Z p . (40) 

The characteristic polynomial of |^0[ ) is r 3 — 2r — 1 = (r + l)(r 2 — r — 1) with a root —1 £ Z p . 
Thus Corollary ["T9l yields a sc-factorization 

t n +l = a 0^n + 0'\tn— 1) Xn+i = —X n + tn+l- 
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where t\ = x$ + x\, ti = x\ + X2- The constant coefficients are a' = a'i = 1 which yield the 
factor equation t n+ \ = t n + t n -\, i.e., the Fibonacci recurrence [26\) . Its characteristic polynomial 
is r 2 — r — 1. Corollary [7PI may be applied as long as IW\) has a unitary solution in 7L V , i.e., a 
solution that never visits zero; Z p contains such a zero-avoiding solution of 126\) for primes of type 
p = 0, 1, 4(mod 5) and infinitely many primes of type p = 2, 3(mod 5); see Jdjj and Let p be such 
a prime and {u n } a zero-avoiding solution of 126\) in 7L V . Repeating the calculations in Example\13\ 
but replacing F n with u n gives the following sc-factorization for |^0| ) 

tn+1 = tn + Sn+1, t± = X\ + Xq 

U n 

U n -1 U2 

Sn+1 = S n , S2 = t2 1\. 

U n Ui 

7 A note on the Poincare-Perron Theorem 

The fact that eigensequences are ratio sequences of unitary solutions recalls the celebrated theorem 
of Poincare and Perron; see [13] and [15], or Section 8.2 of [5]. In the language of eigensequences, 
the theorem may be stated as follows: 



"Let R = C and assume that the coefficients ai^ n in 120\) converge to constants a% as n — )• 
oo ; i.e., i20\) is a Poincare difference equation. Then each eigenvalue of \22}) is a limit of an 
eigensequence of \2(ty . " 



Example 22 Consider the following Poincare equation in the field R of real numbers 

Xn+l = ~X n + X n -x (41) 

n 

The limiting autonomous equation for this is y n +\ = y n -i whose eigenvalues are ±1, i.e., the 
roots of r 2 — 1 = 0. The characteristic equation of m\ ) is 

11 

r„+i = - + — • 42 
n r n 

It is readily verified by induction that the solution of \J^2^ with r\ = 1 may be expressed as 

2n 

T2n-l = 1, r 2n = ~ -■ 

In — 1 

Thus lim n ^oo r n = 1, as expected. Having the above eigensequence explicitly is actually much 
more significant than this convergence result; it yields the general solution of |^i[ ). First, we obtain 
a semiconjugate factorization of |^ 1\ ) with factor equation 

In- 1 

t2n — —t2n-l, *2n+l — — — ^2n- 
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By straightforward iteration 



(2n)l 

t2n+1 ~ 4™(n!) 2 ' t2n+2 ~ ~ t2n+1 ' 

Finally, using the cofactor equation x n +i = r n+ \x n + t n+ i a formula for the general solution of 
(JJ\ ) may be obtained if desired. 

Autonomous difference equations are trivially of Poincare type and each is its own limiting 
equation. Such equations are good for illustrating some aspects of the Poincare-Perron Theorem. 
For instance, to see that not every eigensequence of a Poincare difference equation converges to an 
eigenvalue of the limiting equation, recall from earlier discussion that Equation (|25p has complex 
eigenvalues. Therefore, its real eigensequences cannot converge to such eigenvalues. 

Whether some eigensequences of a Poincare difference equation in topological rings more general 
than R or C (e.g., Banach algebras) converge to eigenvalues of the limiting equation is an interesting 
problem for future discussion. 

8 Zero-avoiding solutions and the ring of quotients 

If R is a commutative ring with identity and no zero divisors (i.e., R is an integral domain) then its 
complete ring of quotients is a field in which R is embedded (see, e.g., |8], Chapter 3). We denote 
this field of quotients by Qr which contains an isomorphic copy of R. The unit group of Qr is the 
set <5i?\{0} of all nonzero elements of Qr which contains i?\{0}. If Qr is not isomorphic to R then 
Qr has an abundance of units that do not exist in R. 

Let us call a sequence {x n } zero-avoiding if x n ^ for all n. In particular, if R is a field then a 
sequence is zero-avoiding if and only if it is unitary. The next result is a consequence of Theorems 
[9] and [18] that reduces the search for eigensequences to a search for zero-avoiding solutions. 

Theorem 23 Let R be an integral domain with field of quotients Qr and assume that the param- 
eters a,j n , b n in (OP are in R for all n and all j = 0, 1, . . . , k. 

(a) If {x n } is a zero-avoiding solution of l2(K) in R then {xnX'^} is a (unitary) eigensequence 
of$MWinQ R . 

(b) If \20\) has a zero-avoiding solution {x n } in R then ^ has a semiconjugate factorization 
in Qr consisting of the pair of equations I116\) and (Jl) with parameters a n = x n x^_ x and a' m n in 
Qr for all m, n. 



The earlier discussion of the Fibonacci recurrence (J26J) illustrates the use of this theorem in a 
familiar case where the integral domain is Z with field of quotients Q. 
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9 SC factorization without unitary solutions 

We now consider a difference equation whose characteristic equation (|23p always has a constant 
solution (eigenvalue) in an arbitrary ring. This difference equation shows that the hypotheses in 
Theorems [9l and 1181 are not necessary for the existence of a semiconjugate factorization. 

Theorem 24 Let R be an arbitrary nonzero ring, let a,b be nonzero elements in R and let {c n } 
be an arbitrary sequence in R. Consider the second-order difference equation 

x n+ i = (a + b)x n - abx n -i + c n . (43) 

(a) The characteristic equation t23\) has a constant solution r = b and |j5p has a semiconjugate 
factorization 

t n +i = atn + c n , x n+ i = bx n + t n+ i. (44) 

(b) If Cn = for all n, i.e., ffify is homogeneous, then every solution {x n } of is given by 
the following where t\ = x\ — bx^: 

Xn = b n x + ^ + b n - 1 + ^"V" 1 j h. (45) 

Proof, (a) In this case the characteristic equation (123|) is 

r n+ ir n - (a + b)r n + a6r n _i = 

The constant solutions of this equation satisfy the polynomial equation (|24p . which is, in this 
case, r 2 — (a + b)r + ab = 0. This evidently has a solution r = b (if R is commutative then r = a 
is also a solution) . Theorems [9] and [18] are not applicable in the absence of an identity and unitary 
solutions, but a sc-factorization of (|43p is readily obtained by rearranging the terms of the equation 
as x n +i — bx n = a(x n - bx n -i) + c n and defining a new variable t n = x n - bx n -\ to obtain the 
desired pair of equations (|44p . 

(b) To establish (j45j) . let c n = and use the first equation in (j44|) to obtain t n = a n ~ l t\ with 
t\ = x\ — bxo. Now iterate the equation x n+ \ = bx n + a n ti and use x\ = bxo + 1\ to obtain 



= b n x Q + (a n " 1 + b n - 1 + ^2b n 



which is (SSI). ■ 
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Example 25 (Boolean rings) Let R be a nonzero Boolean ring defined by the relation r 2 = r for 
all r € R. It is straightforward to show that R is commutative and 2r = r + r = for all r ( consider 
(a + b) 2 which must equal to a + b). A concrete example of a Boolean ring is the collection of all 
finite subsets ofTL (including the empty set) under the operations 

A + B = (A\B) U (B\A), AB = AC\B 

for all finite A,B C Z. This is a commutative ring in which the empty set is the zero element and 
every nonempty set is a zero divisor. 

Theorem ^24\ applies in Boolean rings. From |^5[ j, if c n = for all n then 

x n = bxo + [a + b + (n — 2)6a](xi — bxo) 

bxo + (a + b){x\ — bxo), if n is even 
bxo + (a + b + ab)(x± — bxo), if n is odd 

(a + b)x\ + abxo, if n is even 
(a + b)x\ + abx\ , if n is odd 

where in the last step we used the fact that —r = r for all r. 

The next result shows how formula f)45[) takes a more familiar form in the field of quotients of 
an integral domain. 

Corollary 26 Let R be an integral domain and let a,b G R. Then equation has a semiconjugate 
factorization in R given by (44^ and if c n = then its solution |75[ ) in Qr simplifies to the classical 
formula: 

c\a n + C2& n if a / b 
[nx± — (n — l)6xo]6 n_1 if a = b 

where, using the reciprocal notation 1/u to denote the multiplicative inverse u^ 1 in Qr, 

x\ — bxo axo — x\ 

a — b a — b 

Proof. The first assertion about a sc-factorization is an immediate consequence of Theorem | 
Formula f|45j) is defined in R but it cannot be simplified there. It does simplify in Qr with the aid 
of the geometric sum formula as follows: 

Xn = b n x + b n - x (x x - bx ) (defining a = 1, b° = 1) 

i=i 

b n x + b n ~ l (xi - bx )[(a/b) n - l][(o/6) - l]" 1 if a ^ b 
b n xo + b n ~ 1 n(xi — bxo) if a = b 

Cl a n + c 2 b n if a / b 
[nx\ — (n — l)bxo\b n ~ l if a = b 

where c\ , c 2 are defined as in the statement of the corollary. ■ 



I. ( 46 ) 



20 



Remark 27 Note that in the case a = b fl^ffi ) is defined in R itself. If R is afield then formula fl^ffi ) 
holds in R which is isomorphic to Qr. Care must be exercised in using this formula outside the 
context of fields. For instance, Example ESI shows that (Jffy does not hold in a Boolean ring ( there 
is a mixed product term ab). Further, in rings of functions that are discussed in the next section, 
a, b are parametrized quantities and thus they, or a — b may fail to be units for some parameter 
values. 

10 SC-Factorization in rings of functions 

Difference equations in rings of functions often appear in applied mathematics. Well-known special 
functions such as Bessel functions satisfy recurrence relations that are examples of difference equa- 
tions on rings of real or complex- valued functions (see Examples 1301 and 1531 below). For functions 
from a nonempty set S into a nonzero ring 1Z the operations of addition and multiplication are 
defined pointwise, i.e., for each s £ S 

(/ + g)(s) = f(s) + g(s), (fg)(s) = f(s)g(s). 

Other types of ring operations are possible for functions but we consider only the above pointwise 
operations. With these operations, the set 1Z S of all functions from S into 1Z is a function ring and 
each subring R(S) of 1Z S is a ring of 7£-valued functions on S. Note that R(S) is commutative if 
1Z is. If R(S) contains all the constant functions on S then we usually think of these functions as 
elements of 1Z and thus, think of 1Z as a subring of R(S). In this section we assume that R(S) 
contains all the constants. 

A ring of functions R(S) of the above type is also a function algebra; see, e.g., [8] or [12] . An 
element u is a unit in R{S) if and only if u(s) ^ for all s £ 5. In this case, the inverse of u is just 
its reciprocal 1/u. Since R(S) is closed under addition and multiplication, if the parameters and 
initial values ctj n , b n , Xj : S — >■ M are in R(S) for all j = 0,1, ... k and all n then the solution {x n } 
of ((TJ) is also contained in R(S). 

In the familiar ring C[0, 1] of all continuous, real-valued functions on the interval [0,1] the units 
are functions that are always either positive or negative and a zero divisor is a function whose set 
of zeros has a nonempty interior in [0,1]. The ring of polynomials J~[x] with coefficients in a given 
field J- is a familiar ring that may be viewed as a ring of functions on T (or some subset of it) by 
interpreting the indeterminate as a variable; see [12j . Chapter 4. In particular, if S = [0, 1] and 
T = R then by the Weierstrass approximation theorem T\x\ is a dense subring of C[0, 1] in the 
uniform topology. Various rings of differentiable functions fall in-between T\x\ and the continuous 
functions and share the aforementioned properties of the continuous functions. But larger rings 
such as bounded functions or integrable functions have different properties. 
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Corollary 28 Let R(S) be a ring of real-valued functions on a nonempty set S. Assume that 
a-j,n(s) > for all s £ S, j = 0, 1, . . . , k and all n. If 

k 

Y j a j , n (s)>0 (47) 

3=0 

for all s G S and all n then the homogeneous part of the difference equation |7P has unitary solutions 
in R(S). Therefore, £7]) has a sc-factorization in R(S) that is given by 1116)) and J-?7[ ). 



Proof. Let dj^ n (s) > for all s E S and all n. Choose constant initial values Uj = 1 for j = 
0, 1, . . . , k in (|20l) . i.e., the homogeneous part of dl]). By (13711 . Ylj=o a j,n( s ) > so 

k 

Uk+i{s) = ^2a j:k (s) > 

3=0 

for all s G S. Thus Uk+i(s) is a unit in i?(5") and 

k k 

Uk+2(s) = '^2a j ^+i(s)u k+ i^ j (s) = a , k +i(s)u k+1 (s) + ^0^+1(5) 
3=0 j=i 

for all s G S. If X^=i a i,fe+i( s ) = f° r some s then by (|47|) ao,A:+i(s) / 0. It follows that Uk+2 is 
also positive on S, hence a unit in R(S). Proceeding in this fashion, it follows that u n (s) > for 
all s G S and all n. Thus {u n (s)} is a unitary solution of (120p . By Theorem 1181 the ratios sequence 
{u n (s)/u n _i(s)} is a unitary eigensequence in R(S) so Theorem [9] yields a sc-factorization for ([1]). 
■ 

Remark 29 Nonunitary solutions for \20\) exist under the hypotheses of Corollary [23 because an 
initial function may not be a unit. Further, none of the parameters aj tn (s),b n (s) in Corollary [28\ 
may be units. For instance, the corollary applies to the following difference equation in C[0, 1] 

x n +i{r) = a (1 — sin nirr) x n (r) + br n (l — r)x n -i(r), a,b > 0, r€ [0,1], n>l 

in which ao,n(^) = a (1 — sinn7rr) ; ai >n (r) = br n {l — r) and b n {r) = are nonunits, but for all r,n 

ao,n(?") + ai,n( r ) = a (1 — sinn7rr) + 6r n (l — r) > 0. 

Example 30 (Modified Bessel functions) The second order, linear difference equation 

2n 

x n+ i(s) = — x n {s) + x n -i(», sG(0,oo) (48) 
s 
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is the recurrence relation for the modified Bessel functions K n (s) of the second kind, so-called 
because they are solutions of the second-order linear differential equation known as Bessel' s modified 
differential equation (see, e.g., JJSj/j. In fact, the sequence of functions {K n (s)} is a particular 
solution of ft48\ ) from specified initial values Kq(s),Ki(s). According to Corollaru \28\ a unitary 
solution {u n (s)} of fl^#| ) is generated by any pair of positive functions; e.g., uq(s) = ui(s) = 1. The 
first few terms are 

2 8 4 48 24 2 

U2(s) = - + 1, U 3 (s) = -r H hi, llAs) = ^r + ^-H hi 

S S z S S' 3 s z s 

Now the ratios u n {s) / u n -\{s) define an eigensequence for f^g| ) and yield the sc- factorization 

t n+ i(s) = -— t n (s) x n+1 [s) = -— x n (s) + t n+1 (s) 

u n (s) u n (s) 

with t\(s) = xi(s) — [ui(s)/uq(s)]x (s) = xi(s) — xq(s). Iteration of the factor equation yields 
t n (s) = (—l) n ~ 1 ti(s) /u n -i(s) ; inserting this into the cofactor, summation yields a formula for the 
general solution of fl^ffp in terms of the unitary solution {u n (s)} as follows: 

n-l , v 

/ \ / \ / \ ^ ^ \ ^ J / \ 

X n (s) = U n (s)Xl(s) + } T\ti(s) 



U n (s) 



n-l , 



x (s) + t 1 (s)Y^ ' '' 



^ Ui(s)Ui-l(s) 



Different values of positive functions uq(s),ui(s) yield different formulas but of course, the same 
quantity x n (s). 

In the next result, the coefficients are sequences of constant functions; however, they may exist 
in any nontrivial field, not just the real numbers. 

Corollary 31 Let J 7 be a nonzero field, S a nonempty set and J~(S) a ring of functions from S 
into J- that contains all the constants in T . Consider the following difference equation in ^(S) 

x n +l(s) = a ,nx n (s) + ai )n x n -i(s) H h a k)n x n -k{s) + b n {s) (49) 

where the coefficients aj^ n are constants, i.e., a^ n G J- for j = 0, 1, . . . k and b n 6 ^(S) for every n. 
If the homogeneous part of |^9| ) has a zero-avoiding solution in J 7 (as a subset of F{S)) then |^P| ) 
has a semiconjugate factorization in J-(S). 

Proof. Let yo,yi be constant functions, which may be thought of as elements of J- such that the 
solution {y„} of the homogeneous part of (|49l) that they generate is zero-avoiding. Then y n ^ 
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for all n so that the ratios a n = Un/Un-i are well-defined units in J-(S). Thus by Theorem 1181 we 
have a unitary eigensequence and Theorem [9] yields a sc- factorization for ([49 p . ■ 

Going in a different direction than Corollary [3T] we use Theorem 1241 to obtain semiconjugate fac- 
torizations and solutions of certain second-order difference equations in R(S) to which the preceding 
results may not apply. 

Corollary 32 Let R(S) denote a ring of complex-valued functions on a nonempty set S and let 
fi9,u) n G R(S) for all n > 1. The linear difference equation 

x n+ i(s) = f(s)x n (s) + g(s)x n -i(s) + w n (s) (50) 

has the semiconjugate factorization 

i \ f( s ) + M s ) / \ i \ I ^ — t \ / \ 

t n+1 {s) = t n (s) + w n (s), x n+ i(s) = x n (s) + t n+1 (s) 

in R(S) where h{s) = [f 2 {s) + Ag(s)] l l 2 . 

Proof. Let a{s) = [f{s) + h(s)]/2 and 6(a) = [/(s) - h(s)]/2. Then /(a) = o(s) + 6(a) and 
(/(s) = — a(a)6(a) so (|50p can be written in the form (|43p and Theorem 1241 completes the proof. ■ 
We now apply the ideas in Corollary 1321 and Corollary 1261 to obtain general solutions of linear, 
second-order difference equations in rings of functions. Assume that /(a) and g(s) are real-valued 
functions in the ring R(S) of Corollary [32] and consider the homogeneous difference equation 

x n+ i(s) = f(s)x n (s) + g(s)x n -i(s). (51) 

Let 

S + = { s e S : h 2 (s) > 0}, S- = {s G S : h 2 (s) < 0}, S = {s G S : h 2 (s) = 0}. 

S = S + U S- U So is a disjoint union of the above three sets, at least one of which is nonempty. 
For a G S+ U So the formulas in Corollary [26] yield correct solutions within R(S) to the system of 
equations in Corollary [32] 

If S_ is nonempty and s G S_ then g(a) < and 

h(s) = i^-p( S )-4g(s), a( S ) = ±[f( S ) + i^/-f 2 (s)-4g( S )], 6(a) = <i) 

where the definition of a(s) comes from the proof of Corollary 1321 Switching to polar coordinates, 
the modulus \a(s)\ = |6(a)| = y / —g(s) and 

a ( s ) = \/~ 9( s )[ cos 9( s ) + ism6(s)], 6(a) = y /— 5( s )[cos 6(s) — isin^(a)] 
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where 



9(s) = cos" 



-i m 



2^g{s) 

Repeating the calculations in the proof of Corollary [26] applied to the sc-factorization in Corol- 
lary [32] readily yields 



xi(s) - a(s)x (s) a(s)x (s) - xi(s) 

^n(s) = — == — a (s) + 



a(s) 



a(s) — a(s) a(s) — a(s) 

Dropping the parentheses to reduce notational clutter and changing to polar coordinates, 
straightforward calculation yields 



x n = (-g) 



n/2 



Xq COS n COS 



i / \ , fxo-2x 



+ 



sin ncos 



-i / 



(52) 



Note that if xq,x\ are real valued functions then so is x n for every n. 

Example 33 (Chebyshev polynomials) Consider the linear difference equation 

x n+ i{s) = 2sx n (s) - x n -i(s) (53) 

which is also the recurrence relation for the so-called Chebyshev polynomials T n (s) with the initial 
functions Tq(s) = 1 and I\(s) = s (for each n, T n (s) is a solution of Chebyshev's second-order linear 
differential equation; see, e.a..j23$). If f(s) = 2s and g(s) = —1 then SL = (—1,1), Sq = {—1,1} 
and S + = M\[-l, 1]. For s E (-1, 1), 



h(s) = 2isJ 1 - s 2 , a(s) = s + i\Jl- s 2 

so the formula for the general solution of A53\) is obtained from \52\) as 

, \ I -1 \ sxq(s) — x\(s) . , i \ 

xq{s) cos In cos s) H = sin ncos s) 

\ S 1 * — 1 



x n (s) 

In particular, for n > 1, 



T n (s) = cos (ncos 1 s) 



1< 8 < 1 



a €(-1,1). (54) 



(55) 



which are commonly known to be polynomials. If s G S+ then repeating the calculations in Corollary 
W)[ yields 



a(s) — b(s) 



a(s) — b(s) 



+ Vs 2 -l) + 



xi(s) - (s - yjs 2 - l)xp(s) 
2Vs 2 ^l 

+ (s + Vs 2 -j)x (s)-x 1 (s) ^ _ v /^— [y 
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In particular, 

Tn(s) = \(s+ Vs 2 -l) n + ~ (s - Vs 2 -l) n , s G R\[-l, 1]. (56) 

Finally, if s £ So i.e., s = ±1 then a(s) = b(s) = s so again from Corollary \26\ 

x n {s) = [nxx(s) - (n - l)b(s)x (s)]b n ' 1 (s) 

__ f nxi(l) - (n - l)x (l), if s = 1 
~\ [nx 1 (-l) + (n-l)x (-l)](-l) n , ifs = -l ' 

These relations yield T n (l) = 1 and T n (— 1) = (— l) n /or all n; it follows that formula t56\) also 
holds for s = ±1. Though not immediately apparent, expressions T n defined by ( 155)) and < I56]) taken 
together yield polynomials on the real line. 

11 Summary and future directions 

In this paper we studied semiconjugate factorizations of linear difference equations in rings. For 
the typical (recursive) linear equation this method supplements the standard methods that use 
modules and group actions. 

In some cases, we obtained several sufficient conditions that are not necessary for the existence 
of sc- factorizations. Relaxing or modifying the hypotheses in these cases should lead to broader 
applicability. For example, extending Corollary [28] to include negative coefficients or extending 
Corollary 1321 to include dependence on n will lead, among other things, to results such as Examples 
[301 and [33] on recurrences of many more special functions such as Bessel, Legendre, Hermite, etc. 

Many challenges at different levels of generality remain. A ring with identity may contain no 
unitary solutions of the homogeneous part of (pQ). For example, it is shown in |17j that for certain 
primes, e.g., 2,3,7,23,... the field Z p contains no zero-avoiding (hence unitary) solutions of (i26j) . 
In such cases ([1]) fails to have the linear form symmetry (jlip . However, nonexistence of a particular 
form symmetry is not equivalent to the nonexistence of a sc-factorization; see the remarks following 
the proof of Theorem [9] about the (nonlinear) inversion form symmetry. The question of whether in 
the absence of (11 If) other types of form symmetry exist for which the factor or the cofactor equation 
is linear, remains open. 

If R is not commutative then using different orders of multiplications of coefficients aj ;n with 
the variables x n ^j in ([1]) may result in different equations (with the same set of coefficients). Not 
all the results in this paper extend readily to these variants in the noncommutative cases. 

Finally, linear difference equations may occur in nonrecursive forms; e.g., ao, n x n + oi,n^n-i + 
■ ■ ■ + a k,n x n-k = where the leading coefficient ao, n is not a unit in the ring R for infinitely many 
n and the equation cannot be solved uniquely for x n . Nonrecursive, nonlinear difference equations 
may occur in conjunction with recursive linear equations; indeed, the characteristic equations (|2ip 
and (I23p are nonrecursive polynomial equations. For nonrecursive difference equations even such 
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basic issues as the existence and uniqueness of solutions is not assured. To illustrate, consider the 
linear homogeneous equation 

4x n+ i + 6x n + 2x n _i = (57) 

in the finite ring Zg. Since 4 is a zero divisor in Zg the above equation is not recursive. Hence 
it does not unfold to a map of Zg x Zg in the standard way and semiconjugate factorization is 
not applicable. However, ([57]) does split into a triangular, factor /cofactor pair of equations by a 
simple rearrangement of terms: 4(x„_|_i + x n ) + 2{x n + x n _i) = which by defining the new variable 
t n = x n + x n -i yields 

4t n+1 + 2t n = 0, (58) 
x n+ i = -x n + t n+ i. (59) 

It is evident that every sequence {t n } in {0, 4} N is a solution of the (nonrecursive) first-order 
factor equation (|58p in Zg. The (recursive) cofactor equation (|59[) then generates a unique solution 
of (|57p in Zg for each solution of (158)) : the table below lists the first few terms of a sequence 
{t n } and the corresponding terms of the solution of (f57|) with initial values Xq = 1, x\ = 3 and 
ti = xi + x = 4 



n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 






4 


4 


4 


4 


4 

















4 


4 


4 


4 


4 




Xn 


3 


1 


3 


1 


3 


5 


3 


5 


3 


5 


7 


5 


7 


5 


7 





The same pair of initial values corresponds to infinitely many solutions of (|57p depending on 
the choice of {t n }; e.g., 



n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 




tn 


4 





4 








4 











4 














4 




Xn 


3 


5 


7 


1 


7 


5 


3 


5 


3 


1 


7 


1 


7 


1 


3 





By contrast, the situation is quite different in Zg because 4 is a unit and its reciprocal modulo 9 
is 7. Thus upon dividing by 4, (|5"7]) is reduced to the recursive linear equation x ra+ i+6x n +5x ra _i = 
or equivalently, x n+ \ = 3x n + 4x n _i to which the methods discussed in this paper are applicable 
(though in this simple case the sc-factorization may be quickly obtained via the rearrangement of 
terms mentioned above). These observations indicate that nonrecursive linear difference equations 
and the challenging problems associated with them may turn out to be sources of profound new 
ideas. 
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